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A THEORETICAL STUDY OF THE EFFECT ON EXPANSION- TUBE 
PERFORMANCE OF AREA CHANGES AT PRIMARY AND 
SECONDARY DIAPHRAGM STATIONS 

By Linwood B. Callis, Jr, 

Langley Research Center 

SUMMARY 

An investigation is made to determine the effect on expansion-tube performance and 
operating parameters of the insertion of either a nozzle or an area discontinuity in the 
basic expansion-tube configuration. 

The present analysis is carried out with the area change located at either the pri- 
mary or the secondary diaphragm. During all phases of the facility operating cycle, the 
assumption of a calorically perfect gas is made, thus, analytic solutions are allowed for 
parameters under scrutiny and trends which may be of interest are more clearly 
established. 

From the included analysis, a determination of the performance and operating 
parameters of the device is made. Although certain benefits may be derived from the 
insertion of area changes at the locations described, other considerations make such a 
placement undesirable. 

The advantages as well as disadvantages of the configurations considered are pre- 
sented and discussed briefly. 


INTRODUCTION 

Trimpi (ref. 1) presented a theoretical analysis of the basic expansion tube, delin- 
eating the projected advantages and disadvantages of this device. Of the anticipated draw- 
backs, two prompted the analysis contained in the present paper: viscous effects encoun- 
tered in the acceleration chamber and diaphragm opening problems. 

The first of these problems arises as a result of the long flow lengths at low opera- 
ting pressures encountered in the acceleration chamber. As a result of this combination 
of circumstances, the boundary -layer thickness may be found to be a significant fraction of 
the tube diameter, if this diameter is kept moderately small (4 to 8 inches (10 to 20 cm)). 
This situation may lead to boundary- layer closure at the center line, eliminating any core 
of potential flow, or to large flow gradients resulting from viscous attenuating 



mechanisms. It is, therefore, necessary to find some means to circumvent these effects. 
The most obvious solution to the viscous problem, other than using tubes of large diam- 
eter, is to maintain operating pressures and flow lengths not conducive to severe boundary- 
layer growth. This approach, however, is not practical since, for reasonable test times, 
the expansion tube requires rather long flow lengths. In addition, restricting the range of 
operating pressures negates one of the major advantages of the facility — that is, the abil- 
ity to duplicate flight environments over a wide range of velocities and ambient conditions. 

It is with such considerations in mind that the analysis in the following sections is 
undertaken. The attitude taken is that at some point in the tube, either a nozzle or area 
discontinuity is present which allows the flow to expand into a tube of larger section. 

This area change is located at either the primary (configurations I and II) or secondary 
diaphragm (configuration in). In either case, the acceleration chamber has the benefit of 
an area enlargement, thus, ameliorating the viscous problem. At the same time, the 
diaphragm problem is somewhat diminished in that the use of a nozzle allows smaller 
upstream cross sections for a given test-section size. Thus, at least the primary and 
possibly the secondary diaphragm diameter is reduced depending on the configuration cho- 
sen. The diaphragms being smaller open more quickly and produce more nearly ideal 
wave systems and allow the facility to operate with greater efficiency. 

In review, the nozzle, or area discontinuity concept, is the most direct approach to 
the easement of the viscous and diaphragm problems. However, along with the benefits 
of the nozzle configurations, certain disadvantages must also be accepted: namely, in a 
supersonic flow (which is presupposed behind all traveling shocks) the steady expansion is 
not as effective a converter of thermal to kinetic energy as in a nonsteady expansion 
(ref. 1) and, further, as nonsteady shock waves and/or interfaces pass through the nozzles, 
unsteady adjusting wave systems are necessarily generated which must be reckoned with. 
Such problems and their solutions are investigated and discussed. 

The present analysis and that of reference 2, are complementary and together form 
a rather complete analysis of this type of modification to the expansion tube. From this 
basis, the decision as to which configuration is most feasible can be made and confirmed. 

SYMBOLS 


A cross-sectional area 

A area ratio, A e /Ai 

a sound speed 

l ^ length of driver chamber 
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length of nozzle 

length of intermediate chamber 
length of acceleration chamber 

length of unusable slug test gas in region 2 for configuration DI 
Mach number 

adjusting shock Mach number 

primary shock Mach number 

mass of unusable test gas for configuration HI 

mass of usable test gas in region 5 for configuration in 

pressure 

temperature 

time 

velocity 

adjusting shock velocity 

axial distance measured from primary diaphragm, except where noted 
parameter defined by equation (45) 
ratio of specific heats of test gas, 1.4 

ratio of specific heats of driver and acceleration chamber gas, 1.67 
time interval required for wave type to pass through conical nozzle 
parameter defined by equation (A4) 
parameter defined by equation (39) 
mass density 

ratio of driver pressures, p^ a = 1 0 
mass ratio, m X/ /m5 
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Subscripts: 

A facility with area change 

b,c general subscripts 

e steady-state conditions at nozzle exit 

f flow particle 

i steady- state conditions at nozzle inlet 

m zero length nozzle 

0 constant area expansion tube 

1 configuration I 

II configuration II 

m configuration m 

1 initial conditions in intermediate chamber 

2 conditions behind shock in intermediate chamber 

3 conditions of expanded driver gas 

3e conditions behind adjusting shock 

4 hot driver conditions 

5 test conditions 

10 initial conditions in acceleration chamber 

20 conditions behind shock in acceleration chamber 

+ u + a characteristic 

u - a characteristic 

PERFECT GAS ANALYSIS 

In this analysis, all gases are assumed to be calorically perfect, since this assump- 
tion permits solutions for most of the desired parameters in closed form, which greatly 
facilitates the analysis yet preserves trends of interest. 
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Investigated in the present paper are three modifications of the basic expansion tube 
discussed in reference 1. These modifications are referred to as configurations I, n, 
and IQ. (See fig. 1.) Configuration I (figs. 1(a) and 2(a)) consists of the basic expansion 
tube with a divergent nozzle located immediately downstream of the primary diaphragm. 
The inlet and exit areas of the nozzle are denoted by Aj and A e , respectively, resulting 
in an area ratio A = A e /Ai. The nozzle length is denoted by l n . 

Configuration I is operated as an "ideally expanded" nozzle in that the pressure and 
velocity of the driver gas at the nozzle exit (conditions e, fig. 1(a)) match identically 
those in the shock initiated flow (region 2). Thus, no waves are necessary to adjust the 
nozzle and shock flow to a common interface condition. Pressure ratios are taken to be 
such that a value of = 1 is insured. 

Configuration n is physically identical to configuration I, the difference being in the 
mode of operation. Configuration n is operated in the over-expanded mode which is to 
say p e < P 2 , u e > U 2 . Hence, an upstream-facing shock wave is required in order to 
bring conditions e to state 3e, which is identical in u and p to state 2. (See 
fig. 1(b).) Forward-facing adjusting waves in this system are not permitted. 

For configuration III (figs. 1(c) and 2(b)) the nozzle is located immediately down- 
stream of the secondary diaphragm. This location results in the primary as well as the 
secondary diaphragms being reduced in size and the acceleration chamber being increased 
in diameter. The anticipated wave schematic for this configuration is shown in figure 1(c), 
which, along with figures 1(a) and 1(b), may be compared with the wave schematic of the 
basic expansion tube (fig. 3). Some experimental and theoretical work on this configura- 
tion has been done at Arnold Engineering Development Center. 


Driver Requirements for Configuration I 

One means of analyzing the present configurations is by comparing their operating 
parameters to those of the basic expansion tube for identical test states. Before begin- 
ning, however, the following relationships which are basic to steady and unsteady flows 
and used repeatedly are presented along with their regions of applicability: 
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(Riemann, invariant for unsteady (la) 

isentropic constant area flow) 

(steady homenthalpic flow) (lb) 


(isentropic) (2a) 
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(isentropic) ( 2 b) 
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where the subscripts b and c refer to any two states along an isentrope. 

Consider now the driver pressure required for the basic expansion tube for a given 
test-section Mach number The required relation is 
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which is equivalent to equation (36) of reference 2 for A = 1 and = Mg. In equa- 
tion (3), the strong shock limit for M 2 , 


m 2 = 


y(y - 1) 


(4) 


has been used. Henceforth, it is assumed that all shock waves appearing in this analysis, 
with the exception of the adjusting shock in configuration II, are strong shocks. Equa- 
tion (3) is shown plotted in figure 4 for (a^/ag) = 5, 10, and 20 and for ~ = 1.67 for a 
helium driver. 

For configuration I, the same parameter P 4 /P 5 may be expressed in terms of 
ratios as follows: 




(5) 


These ratios (eq. (5)) may be evaluated by assuming that conditions at nozzle inlet are 
sonic and by using equations (la), (lb), and (2a) with the interface relation P 2 = p . When 
these values are substituted into equation (5), the result is 
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( 6 ) 



Implicit in equation (6) is the assumption that waves generated by the starting process of 
the nozzle are negligible. This assumption is discussed in a later section. 

At this point, having determined p 4 so that P2 = p e , the sound speed a 4 neces- 
sary to insure that u e = u 2 must be determined. Let 

u 2 = u e = M e a e (7) 

where 


o _* a e_ u 2 
e 4 a 4 M e 


( 8 ) 


with 


u 2 



0) 


as determined from equation (la). Substituting equation (9) and a e/ /a 4 , determined from 
equations (la) and (lb), into equation (8) and simplifying gives the required condition for 

(^/^j as 
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M e \|r(r - 1) 


( + £-l_±Mg 


r - 1 • 
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( 10 ) 


This relation is a function only of 5 and Mg and is plotted in figure 5. As is apparent 
from the figures, the required value of (a^^ag) increases with Mg. More interesting, 
however, is the fact that the required sound speed ratio drops sharply as A increases to 
a value of 10. Above A = 10, the driver sound speed becomes relatively insensitive to 
changes in A. In fact, as A (i.e., M e ) becomes larger, (a 4 /a 5 ) approaches a limit 
independent of A, given by 
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( 11 ) 




Also indicated in figure 5 are the values of a 4 /ag at which T 4 /T 5 = 50, 25, and 10, the 
gases being calorically perfect helium and air in regions 4 and 5, respectively. It 
becomes apparent then that, for the helium air system with A = 10, T = 4200° K for 
T5 = 300° K. These, with present arc discharge techniques, are relatively modest 
requirements. 


With the use of equations (3), (6), and (10), the parameter <f>j given by 


(V H 
(VU 



Jr 

7-1 


( 12 ) 


is determined and is indicative of the relative change in the driver pressure required as 
a result of the nozzle addition to the basic expansion-tube configuration. Equation (12) is 
shown in figure 6 and is a function only of A, this being a result of using the same value 
of a 4 /a5 in the expansion tube as in configuration I. As seen from figure 6, configura- 
tion I suffers large pressure disadvantages for nozzles of large A. These relative dis- 
advantages become more severe if values of a 4 /a.g, larger than required by configura- 
tion I, are used in the expansion tube. 

Several comments should be made at this point. First, in figure 6, for 
A = 1.0(log 10 A = 9) the basic expansion-tube process is not represented in general, even 
though the cross-sectional area is invariant. The process indicated in this fashion is the 
limit 


M e - = 1.0 (13) 

representing expans ion-tube operation with the restriction that the driver gas is expanded 
only to sonic conditions. This is not a general restriction on expansion-tube operation 
and henceforth, for all configurations, curves labeled A = 1.0 represent the limit 
described by equation (13), whereas those labeled A = 1.0 denote the general expansion- 
tube cycle. Any exceptions to this are noted. 

The asymptote indicated by the vertical line noted in figure 4 represents the expan- 
sion of the driver to the limiting state. In the limiting state, pg/p 4 — 0 and p 4 
moves to infinity in an effort to raise the level of Pg to the nonzero value of pg and to 
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maintain interface requirements. The net result for finite Pg is to drive ^> 4 /p 5 ^ — «>'• 
The limiting states for (a 4 /a 5 ^ = 10 and 20 occur at values of Mg given by 


Mg = 


y - i\ a 5 




(14) 


and are well out of range for the conditions of interest herein. 

Driver Requirements for Configuration II 

From the preceding section, it becomes clear that configuration I will generally 
operate at a disadvantage, relative to the expansion tube, with regard to driver pressure. 
For A = 10^, the expansion tube is capable of operating at pressures 4 to 5 orders of 
magnitude lower than that required in configuration I. Driver sound speed requirements, 
however, are modest, and at low values of Mg, quite low. However, if ( a 4 /a g)_ is 
raised beyond those values required for a perfectly expanding nozzle, one would expect a 
flow requiring a backward-facing shock (fig. 1 (b)) and possibly other unsteady waves in 
order to establish the interface conditions between 2 and 3e (fig. 1(b)). This adjusting 
wave system would imply driver pressures lower than those required for configuration I. 

Consider first the relative magnitude of the nondimensional velocities Ug/ag and 
u e /a 5 and of the pressures p 2 /p 5 and P e /P 5 . From equations (la) and (4) applied 
across regions 2 to 5, the following expression is obtained: 


u 1 + 

u o 2 0 
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y 


2y 


- 1 Vy(y - 1) 


(15) 


The quantity u e /ag may be written as 


u c 


5/j\ 4 


a iV a eV u e 


(16) 


When substituting equations (la) and (lb) into equation (16), it may be written as 


u p /a 4 \ 

_ - _ M 
a. “ lap- L AVA e 
5 \ 5/ a 


y + 1 


i + ^1 m, 


- 1/2 


(17) 
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Results of equations (15) and (17) are shown graphically in figure 7 in the form 

. The ticks on the curves indicate the values of (a^/ag)- at which the 


u e 

= f 
a 5 


4> a 


nozzle would be perfectly expanded (iig/ag = u 2 / a 5 ) f° r that particular value of A and 

Mg. If for a given value of A and Mg, (a^yagj- exceeds this value, the nozzle will 

be overexpanded since u e /ag > u 2 /ag. From figure 7, it is apparent that for 
(a^/ag)^ 2 11 , a modest sound speed requirement, the nozzle will be overexpanded for 

M 5 = 50, A S 10. 

When applying equations (la) and (2a) across regions 2 and 5, the nondimensional 
pressure p^/Pg is 


/l + 


^“5' 


Jr 
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V. 1 + l^ 

With the aid of equations (la), (lb), and (2a), the quantity Pg/Pg is given by 


y - 1 


(18) 


V p 5A p iA p 4y 


P 4 
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(19) 


"Mg 2 ' 


The result of the combination of equations (18) and (19) is shown in figure 8 where the 
l- p 4\ 

function P e /Pg = f|A,— j i s s h° wn plotted. Superposed on this plot are a family of verti- 


cal lines ? each representing a value of Mg. The intersection of the vertical lines with the 
lines of constant A define the values of P 4 /P 5 required for perfect expansion for that 
particular value of A and Mg. Values of P^/Pg lower than the perfect expansion value 
produce an overexpanded flow requiring a rearward-facing shock to adjust Pg/Pg to the 
level of p 2 /Pg. 


Figure 8 also makes clear the fact that, for the higher values of Mg and A, values 
of p^ required to produce a perfectly expanded flow are too high. Operation in the over- 
expanded mode appears to be necessary. 


The question to be considered is what type of adjusting waves is required. At this 
point, lengths and pressures involved are assumed to be such that the flow between the 


10 


primary shock and interface may be considered uniform. Driver conditions are then 
selected such that an upstream-facing shock wave, moving into the expanded driver gas, 
produces the adjustment required for the overexpansion. This process is indicated in fig- 
ure 1(b) where 3e denotes interface conditions in the driver gas. 

Since p^ = Pg g , the parameter (P4/P5) may be written as 


l p J\ ( P M- lYM^ 

\ P 5/ n ~ \ P i/v P eA P 3e/\ P 5/ 


By dividing by , the parameter 


*n is defined as 


( 9 *l P s)g 

W) 0 


( 20 ) 


(21) 


The ratio Pg e /P e in equation (20) is given from normal shock relations as 

(22) 

where 


M 


es 


u e + u es 


(23) 


The adjusting shock velocity u es relative to the tube is assumed positive in an upstream 
direction. 


Appendix A shows that the required value of M es 
the expression 


for such a situation is given by 
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where 


u e/ a 5 - u 2/ a 5 
(a e /a 4 )(a 4 /a 5 )_ 


and depends on 



Mg, and A. 
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In figure 9, the shock Mach number of the required adjusting shock is shown as a 
function of ( a 4 / a 5 ^> Mg, and A. This plot was prepared with the use of equations (la), 

(lb), (2a), (15), and (16) to evaluate £. Equation (A5) was then used to determine M eg . 

As seen from figure 9, the adjusting shock Mach number decreases with increasing 
Mg. In fact, for (a 4 /ag^_ = 5 and 10, M eg assumes the value of unity thus indicating 
a perfect expansion within the Mg range considered herein. The value of Mg at which 
this occurs is given by the inverse of equation (10). For values of > 10, Mg is 

out of range. 


With the aid of equations (1), (2), (3), and (22), <E>jj may be determined for a given 
value of Mg, A, and ( a 4 / a 5 ^* Values of this parameter are presented in figure 10 for 

the several values of (a 4 /agj_ and A. 

The conclusion that significant reductions in driver pressure may be achieved for 
A £ 10 2 as a result of the overexpanded configuration may be drawn from figure 13. The 
costs of such gains are the larger values of ^a 4 /agj_ required (when compared with con- 
figuration I). These values lie within the expansion-tube range, however, and are 
tolerable. 


One final point of interest is noted. 


For the curves 


for ( a 4/ a 5) 5 


= 5 and 10, as 


those values of Mg are approached at which M eg — 1.0 (fig. 9) then configuration II 
becomes perfectly expanded. This condition is reflected in figure 10 by the upswing of 
those particular curves terminating at the proper value of Mg given by the inverse of 
equation (10). The points of termination are joined by the near vertical curves shown in 
figure 10 labeled "Perfect expansion limit." 


There must now be answered the question of when the adjusting shock wave is 
trapped in the nozzle. Such a situation would result in an unstarted nozzle, possibly gen- 
erating strong upstream waves, and is not permissible. From equation (23) u eg , the 
wave velocity relative to the facility, is given by 


u „ = M„a„ - u 
es es e e 


(24) 


Since u is taken as positive in an upstream direction, it is required that u „ g o. 

co GS 

This inequality insures that the nozzle flow is not disturbed by the adjusting wave. Equa- 
tion (24) may then be written 

u es a e - u e < 0 ( 25 > 


12 


I 



or 


M es < M e (26) 

Equation (26) is now the requirement to be satisfied. From equations (1), (2), (15), (16), 
(A4), (A5), and (26), values of Mg required to avoid trapped shocks were determined for 
a range of A and ^a^ya^— . The results, shown in figure 11, ideally indicate the lowest 

a 5 )^. Applica- 
tion of these lower limits is indicated in figure 10 by the lines, labeled "Trapped shock 
limit," terminating the curves at the lower values of Mg. Before proceeding, it 

should be pointed out that the possibility of a third limiting curve exists. This limit (not 
shown here) would in practice occur at conditions where the adjusting shock causes 
boundary- layer separation possibly generating downstream waves that render the flow 
useless. It is believed that this limit occurs at values of Mg larger than the trapped 
shock values; thus, the already small operating regime is further truncated. 


values of Mg which may be tolerated for 




Driver Requirements for Configuration in 

Configuration HI, investigated in a preliminary manner at the Arnold Engineering 
Development Center, whose wave and physical schematics are shown in figures 1(c) and 
2(b), is assumed to operate by first expanding the shock-processed gas 2 steadily through 
the inserted nozzle to conditions e then unsteadily by the u - a expansion fan to test 
conditions 5, with interactions due to starting phenomena being neglected. The following 
analysis proceeds from this assumption. 

The driver pressure parameter ^p^^/p^j^ may be written in terms of the ratios 


since 




p 3 = p 2 = Pi 


(27) 


(28) 


By utilizing the interface conditions along with equations (la), (lb), (2a), (4), and (28), 
equation (27) may be expressed as 
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( 29 ) 


Dividing equation (29) by equation (3) ^with ( a 4/ a 5)^ = ( a 4/ a 5 ) Q j resu lts i 


in the ratio 


m 


<t> 


m 


( p 4/ p 5) m 

( P 4/ P 5)„ 



(30) 


This ratio is shown graphically in figure 12 for various values of the independent varia- 
bles. As is evident from this figure, there are several limitations, denoted by the labeled 
vertical lines, which restrict the range of operation of this configuration. The indicated 
asymptotes represent limit lines to which the curves for A = 10, 10^, and 10^ for 
( a 4/ag)- = 5, tend. This limit is a result of the limiting state being reached by configura- 
tion III resulting in P4/P3 — °°. A comparison of figure 12 with figure 4 indicates that 
this limiting state is achieved at a smaller value of Mg for configuration HI than that for 
the basic expansion tube. The last limitation, indicated in figure 12, represents the Mach 
number to which the inlet flow expands when passed through a nozzle with A = 10^; con- 
sequently, this is the lowest value of Mg which can be achieved with this configuration. 
This limitation presupposes an inlet Mach number given by strong shock approximations 
(eq. (4)) with y = 1.4. 

Comparison of figure 12 with figures 6 and 10 indicates that in some operating areas, 
configuration EH is the facility with the least penalty in regard to the driver pressure but, 
nonetheless, there remains a penalty of several orders of magnitude in p^ when com- 
pared with the constant-area tube. Further discussion is presented in later sections. 
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Intermediate-Chamber Pressure for Configurations I, n, and HI 


The pressure ratio across the primary shock may be given as 


P 

P 


2 

1 





2 


( 31 ) 


with the strong shock assumptions. If a nozzle is assumed to be present in the configura- 
tion, the quantity P^/Pg is 


Pi 

P5 


Pe P5 


2y 

y + 1 


M c 


(32) 


By applying the analysis in appendix A, it may be determined that 


M Sl 


y + 1 2 

2 )|y(y - 1) 



(33) 


where the strong shock approximation has been employed. The intermediate- chamber 

parameter — (— ) , evaluated with the aid of equations (la), (lb), (2a), and (33), is shown 
p 5 \ 1/ 

in figure 13 as a function of A and and increases with both. 

At this point, it should be noted that once the primary shock wave is set in motion, 
the problem is one of a constant area for configurations I and II. Hence, the curve for 
A = 1.0 in figure 13 is valid for these configurations as well as for the expansion tube. 
However, since the area change does affect configuration III, curves for various values of 
A are shown in figure 13. 


Acceleration-Chamber Pressure for Configurations I, n, and in 


p in /a,- \2 

The acceleration-chamber pressure parameter ma y b e determined for 

p 5 \ a 10/ 

the present configurations by using the analysis of appendix A. The resulting relation 
given by 


/ \ 2 

P 10 a 5 j 2 

p 5 \ 10/ y(y+l)M 5 2 


(34) 
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is shown in figure 14 in inverted form. Vacuum requirements in the accelerating section 
are reduced with decreasing 

Finally, equation (34) is valid for all configurations including the expansion tube 
since once the secondary shock has been set in motion the process is a constant area one. 

Component Lengths 

The determination of the component lengths and their required variation with the 
independent operating parameters to insure maximum operational efficiency are discussed 
in this section. In order to secure explicit solutions for the length ratios of interest, it is 
necessary to neglect nozzle starting effects as well as variation of interface and shock 
velocity in nozzles. These assumptions, though contrary to the physical situation, allow 
solutions from which valuable information, not easily obtained otherwise, can be 
determined. 

Configuration I.- For configuration I the optimum lengths occur when the reflected 
head of the driver expansion, the driver-driven gas interface, and the head of the second- 
ary expansion which processes the test gas coalesce at a given point on the x-t diagram 
(fig. 1). Selecting a driver longer than that necessary to achieve this wave configuration 
results in increasing the energy required in the driver with no increase in testing time. 
Decreasing the driver lengths, however, results in a truncated test time. Hence, the wave 
schematic shown in figure 1(a) represents the most efficient operation of the device and is 
used as a criterion for determining relative component lengths, that is, l^jl 

and Is^jlsy Th e analysis is greatly facilitated with the use of sketch 1 which follows: 


t 
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According to sketch 1, the following equation may be written: 


t 


e 



where, according to reference 3, 


y+1 ) 

tq _ / a 4\2(r _1 ) 

h = w J 


( 35 ) 


(36) 


Since conditions i are sonic, equations (la) and (36) may be combined to give 


y+1 





(37) 


The last term on the right of equation (35) may be given by 


where 






(38) 


(39) 


may be evaluated from reference 2, for a conical nozzle, with A and y known. The 
time tj. is now determined as a function of the known parameters. In order to do this, 
it is necessary to assume that the shock and interface velocities are constant despite the 
presence of the nozzle. This assumption, which is discussed in the section "General Con- 
siderations," allows the following equations to be written: 


Xf = u 2 t f 


, J * Sl +*n o ' 
" *=l + l n + (“2 - a 2 )rf - rfn 


(40) 
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Solving equations (40) for tj with the use of equation (9) gives 


t f - 


Zsj + Zn 

fl + 

2 1 

\(y ~ 1 

a l a 5 

? + V M J 

\7 + 1 


( 41 ) 


Finally, in order to insure coalescence of the desired wave fronts, the following equation 
must be satisfied 


u 2 tf = Z n + (u e + a e )(tf - t e ) 


(42) 


Substituting equations (10), (35), (37), (38), and (41) into equation (42) gives the following 
expression for 


(*o Ni 



(43) 


where 



relation 


is the value of 



for 


r-^- = 0 and is expressed by the functional 
* s i 
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y + 1 (i , y - 1 m 
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y- i, 2 

y(y - 1) 

M e l 

y(y - 1) 

[2 K 2 M e l 


y + 1 y 4- 1 \ 
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y + 1 


(M e + 1; 


y^x\2(r-l) 

k 2 / 


(44) 


The quantity /3 is given by 



with £ — as A 1.0; hence, B — 0. For A = 1.0, equation (44) reduces to 

+ M e + 1 
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! 2 (V + A r - 1 , 2 , |y(y -~l) 

y(y - 1)\ 2 / y + 1 y + 1 i 2 

4 ^ 


and assumes a value =0.431 for y = 1.4 and y=1.67. 

Equation (43) is presented graphically in figure 15 for various values of A and 
Inj^s^ The parameter (idjls^J decreases with increasing A and decreasing 

l n jl s ^. For purposes of comparison, the length ratio (^d/ls^j- is shown super- 
posed in figure 15. This ratio was computed with the aid of the following equation: 


l/M 1 + [2; 


— : m, 

2 5 


y + 1 y + 1 1 


- 1 2 
2 \jy(y - 1) 


( y - 1 
1 + 


r - 1 a 4 


The derivation of equation (47) is not presented because it is similar to that of 
equation (43). Computations using this equation were made at the same value of 
(a 4 /a 5 ) o as required for configuration I. (See eq. (10).) As would be expected, 

(fd/ l s 1 J I < (fd/ZsiJ^ i q ^ or ^ > and s l = In Edition, it is clear from 

figure 15that (idjls^ < (^ d ^s 1 ^_ 1 Q for A > 40 and l n /lsi * 0. 

The length of the acceleration chamber Z s ^ is determined such that maximum test 
time is available. The governing equation may be determined from equation (63) of ref- 
erence 2 and is given by 


IrA l st A y + 1 ' 53 y 


' s i \ l *i 7+1 \i + 2^i M , 


3-r 
2(y- 1) 


where account has been taken of present nozzle length effects and the entropy gradients in 
region 2 have been neglected. Values of M 2 are given by equation (4). Equation (48) 
is shown graphically in figure 16 and is valid for arbitrary values of A for configura- 
tions I and II as well as the expansion tube. As can be seen ~ 0(10 - 103), hence, 
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considerable lengths may be involved for reasonably sized intermediate sections. A more 
complete discussion of this aspect is given in reference 1. 

Configuration II.- The procedure for configuration I is now repeated for configura- 
tion II. The following sketch (sketch 2) is helpful: 



The procedure is the same except that account must now be taken of the fact that the 
reflected head of the expansion encounters the backward facing shock and alters its speed. 

Writing an expression for the axial coordinate of the adjusting shock wave and 
equating it to that for the u + a characteristic results in the following expression: 


* l n + (“e + a e)(‘f s ‘ l e) 


(49) 


Solving equation (49) for the time at which the two waves intersect t, gives 

Ic 


+ ^ n " ( u e + a e)^e 

if - - 

“es + ( u e + a e) 


(50) 


The negative sign is required on the left-hand side of equation (48) since u g is positive 
in an upstream direction. 

If the displacement equation of the interface 

x = Ugt (51) 

and of the u + a wave after intersection with the adjusting shock wave 
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(52) 


x = - u es% + ( u 3e + a 2e)( t " %) 

are written, the conditions for most efficient performance of the device may be deter- 
mined by setting t = tj and equating equations (51) and (52). This operation and the sub- 
stitution of equations (35), (37), (38), (41), and (50) results in 


where 






y - 1 

y + 1 


+ 


2 

y(y - !) 

M es+ 1 

y + A 

2 

M 3e + 1 


K 


+ 1 


y + 1 


y + 1 

W^) 


(53) 


(54) 


Shown plotted in figure 17 is the parameter 



(M e + 1) 


y + 1 


y+1 

Wv\ 



from which (idj^s^j ma y ke determined. Figure 17 includes, for comparison pur- 
poses, a curve representing normal expansion-tube operation, denoted by 5 = 1.0. For 
this case, the ordinate is simply since = °* The quantity 


(idjls j q was determined from equation (47). 


From figure 17, trends become apparent and it is seen that the addition of a nozzle 
to the expansion tube and its operation in the mode of configuration II results in a driver 
of shorter length, this length decreasing with increasing A and decreasing (a 4 /a 5 )_. 

The driver length for configuration I is smaller over most operating ranges than that 
for the expansion tube or configuration II. (Compare figs. 15 and 17.) This result is 
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primarily due to the lower sound speeds required by configuration I coupled with the rela- 
tive efficiencies of the steady and unsteady supersonic expansions in configuration I and 
the expansion tube, respectively. 

It should be pointed out that in figure 17 the curves for A = 10, 10^, and 10^ are 
terminated by a dashed line at the values of Mg which represent the perfectly expanded 
states of configuration I. 

Configuration III .- Consider now configuration HI. From figure 1(c), it is apparent 
that the driver length should be tailored so that the expanded helium driver gas and air 
interface and the head of the reflected nonsteady driver expansion fan coalesce at the 
nozzle entrance. A larger or shorter driver will move the operation of the device off the 
point of maximum efficiency. 

The analysis necessary to determine these lengths is similar to that required for the 
other configurations except that there is no nozzle to contend with and the point of inter- 
section of the interface and expansion wave is fixed. The equation of motion for the inter- 
face is given by equation (51) which is 


x = Ugt 

whereas that for the u + a wave resulting from the reflected expansion is 


x = x 3 + (u 3 + a 3 )(t - t 3 ) 


(55) 


where Xg and 
uniform region 


tg define the point at which the first reflected characteristic enters the 
3. Note that Xg is given by the relation 


x 3 = ( u 3 ' a 3) t 3 (56) 

Setting x = l s in equation (50), the time of interface arrival at the end of the driven sec- 
tion tj is determined. Using this time, equation (56) and x = Z s ^ in equation (55), the 
following expression is obtained: 


z s x “ 2u 3 t 3 


(57) 


Applying equations (36) 
used in equation (57) to 


across the unsteady fan, region 4 to 3, tg may be determined and 


yield the desired expression for 




VIH 


, which is 


22 



( 58 ) 



using the strong shock approximation, steady and unsteady flow equations. 

Hence, with this information, (l d jl s ^ ^ is found as a function of A, Mg, and 
Equation (58) is presented in figure 18 for a range of parameters, the dashed 

line labeled A = 1.0 representing the parameter for normal expansion tube operation. 

The four vertical dashed lines of figure 18(a) represent the values of Mg at which 
(l — 0, a result of achieving the limiting state in the driver expansion. From fig- 

ure 18, it is apparent that significant reductions in the driver to intermediate- chamber 
length ratio may be realized when compared to the expansion tube. The magnitude of 
these benefits decreases with increasing ^yh.gj_. 

It is now desirable to determine the required length of the acceleration chamber, 

The analysis necessary is more easily understood after a consideration of figure 19, 
an x-t representation of the unsteady wave processes occurring in the nozzle and accel- 
eration chamber of configuration III. Shown are the primary and secondary shock waves 
and interfaces, the unsteady u - a expansion fan (resulting in part from the motion of the 
secondary interface through the nozzle), a particle path, and the u + a characteristic 
generated by the entrance of the primary interface in the nozzle. 

Consideration of this figure reveals two major points of interest. The first is that 
particles entering the nonsteady fan prior to its exit from the nozzle may not be considered 
as suitable for testing purposes. Such particles, entering the fan at different times within 
the nozzle, experience different thermodynamic histories and hence, upon exiting from the 
expansion, have variable thermodynamic and flow properties. This flow is denoted by con- 
ditions 5*. Particles entering the unsteady expansion after its exit from the nozzle, when 
allowed to move entirely through the expansion, experience identical flow and thermody- 
namic histories. This mass flow is used for testing purposes and is denoted by 
conditions 5. 
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The second point of interest is simply that for a given value of l s ^, or 

there is a maximum nozzle length, which if exceeded allows no usable test flow at the test 
section (end of This occurs when the nozzle length is such that the first reflected 

expansion characteristic, traveling with a velocity u + a, coalesces with the first usable 
test particle and the last expansion characteristic of the secondary fan at the end of Zgg. 
It is clear that to produce usable test time nozzles shorter than this limiting length must 
be used. The situation for zero test time ( r g = Oj is indicated in the following sketch 
(sketch 3): 



Solutions for this limiting or maximum nozzle length are presented later. 

Now for a given value of l n , ignoring for a moment the leading u + a wave, the 
question to be answered is what value of l is required to produce usable test time. 

As seen from figure 19, if l is excessively short, those particles suitable for test 
purposes may not exit from the fan before passing the observation station. Consequently, 
no uniform region of flow exists. An exact analysis of this problem would require a dif- 
ficult unsteady characteristic solution beyond the scope of this paper. In order to make 
possible analytical solutions, the following approximations are made: 

(1) Conditions 5 ~ 5 f for this aspect of the analysis 

(2) Interactions between the secondary shock and its trailing interface are weak and 

may be neglected 

(3) Effects of the above interactions are not transmitted to regions 5' or 5 

(4) The unsteady expansion for, after exiting from the nozzle, may be treated as 

centered 

These approximations are discussed in a later section. 
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The question raised previously with regard to l s ^ is now approached by utilizing 
a mass-flow analysis. Let the mass of gas contained in the slug of length l x (fig. 19) 
which enters the unsteady expansion wave during passage through the nozzle be unaccept- 
able for test purposes. This unacceptable mass of gas is compared to the mass passing 
the end of the acceleration chamber in the time interval t’ between the arrival of the 
interface and the last u - a characteristic of the secondary expansion. From a compar- 
ison of these masses for specified parameters, that is, A, Mg, and Zsg/Zn, ^ can 
determined whether any useful test time exists and, if so, how much. 

For clarification of the following derivation, consider the following sketch: 



It is first necessary to determine l x . The quantity A may be defined as 



Solving for l x with continuity and the strong shock approximations gives 


(60) 


S *TT U 2 - (*- 




(61) 


where the quantity (A_ - A^ may be written as follows 

A. - A( = - ? f )|a (62) 

with | and ^ being evaluated once again for a conical nozzle with the use of refer- 
ence 2. With equations (61) and (62) the mass of slug is given by 


I 
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(63) 


m x - Pl A i ^TT U 2(«- - *f)i? 


The mass passing the end of the acceleration chamber in the interval r T is 


m 5 = ^5 u 5 A e T ' 


(64) 


where 


H. 


T* = — ^ 


a 5 M 5 (M 5 - 1) 

If i p = m x/ /mg, the following expression may be written 

p 2 P _e 

°e p 5A^s 2 /V 


„ "2 "e 1 1 ,» , \ a 2 | 2~ / M 1N 


/ a e My - i) 

Substituting equations (la), (lb), and (2b) into equation (66) yields 


(65) 


(66) 


^s 


2^ = —I 


1 + 


y 




l n A\r(y “ !) 


y + 1 

r 


V 1 + 


( M 5 - *)(«- - «f) 


(67) 


Equation (67) is presented in figure 20. For most of the operating range, the param- 
eter {ls 2 ^n) x l / is greater than unity by several orders of magnitude. Recalling the 
definition of i//, it is realized that in order to have any usable test time, must be 

a value such that if/ is determined to be less than unity. The conclusion is then, that for 
usable test time, the parameter l s 2/^n must be of the order 10^ to 10^ for 30. 

For realistic nozzle lengths, this results in values of l s ^ which are untenable. From 
figure 20 the observation is that the situation worsens for decreasing A, the curve for 
A = 1.0 representing not the expansion tube but the solution of equation (67) as A — 1.0. 

With regard to the effect of the assumption that state 5 <= 5 f , it is noted from equa- 

y+1 

tion (67) that (Is2^nj4 / ~ (Mg)!'”*'; hence, an overestimation of by 30 percent would 

raise (l S 2 ^n)^ a factor of 5. It is unlikely that conditions 5' would deviate by 
more than 30 percent from 5; however, if such errors were present, l s would still be 
excessively long. 
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In order to make more clear the reason for the large lengths involved, the move- 
ment of a flow particle through a simple centered unsteady expansion is examined. In the 
following analysis, subscripts e and 5 denote, respectively, the conditions of the par- 
ticle when entering and exiting from the expansion represented in the following sketch: 

Particle 

path 


x 

0 x e x 5 

Sketch 5 

The following derivation of the governing equation for the particle motion 





\ 1 + z i 



( 68 ) 


is similar to that of equation (16) of reference 7 which describes the motion of a P 
characteristic through a u - a expansion. Equation (68) is shown in figure 21 for a 
range of parameters. Those values of M e selected represent those existing at the exit 
of a nozzle with area ratio A, y - 1.4, and inlet Mach number, Mj = 1.89. From fig- 
ure 21, this ratio is seen to be extremely large indicating that flow particles moving 
through an unsteady fan require considerable distances before emerging into the uniform 
flow region 5. This is especially true when Mg » 1. It is precisely this effect that 
requires the large values of l s in order to achieve a usable test time in 
configuration III. 

At this point, the evaluation of the effect that the u + a wave system has on the 
length parameter l n [l s ^ is necessary. Briefly, l n jl s ^ is evaluated under the con- 
straint that it be a value such that the zero test time situation discussed previously and 
pictured in sketch 3 prevails at all conditions. Larger values of l n jl s ^ permit no 
usable test time regardless of Zg^. Shorter values allow test time to establish itself. 
Consequently, evaluating l n jl s for a given value of Z s ^ and Tg = 0 establishes a 
useful maximum value for Z n /Z s ^ which must not be exceeded or equaled if finite test 
time is desired. From the derivation in appendix B the relationship expressing this 
limiting value of is given by equation (B7) which is 
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Mr - 1 l 1 + 21 y^ M 5 


y+1 

2(y-l) 


M e — 1 \ n y-1., 
e \i + — M e/ 


- 1 


m 5 - 1 | 


“e- 1 ' 

\ 1 + r l lM e 


2_ 

-1 


;(«- - « + )K 2 


- 1 



and depends on A and Me. The parameters £ and ? are evaluated according to 

/ ^ \ ^ "" "** 

reference 2, whereas (£ 32 /^ nj’/' is determined from figure 20. The tilde notation for 
\p indicates 4/ = 1 . 0 , the value assumed when the first usable particle exists from the 
unsteady fan precisely at the end of 1$^ the situation prescribed herein. 

Equation (B7) is shown in figure 22, with the curve for A = 1.0 representing the 
limit of equation (B7) as M e — Mj, not expansion-tube operation. It is clear from this 
figure that the maximum permissible nozzle lengths are, for the most part, small when 
compared to In fact, is so small that, for reasonable values of Z Sl , 

the steady flow expansion characteristics may often extend past the nozzle exit before 
processing the center- line flow. 

A second consideration is that the flow, if expanded through a large A in such 
short lengths, surely is not one dimensional or even a close approximation to it. Conse- 
quently, a careful evaluation of the effect of the recompression shock waves generated at 
the juncture of the conical nozzle and the accelerating section must be made. These waves 
could easily be detrimental to the test flow. 

DISCUSSION 
General Considerations 

At this point, the major assumptions employed in the foregoing analysis will be 
reiterated and briefly commented on with regard to their purpose and validity as follows: 
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(1) All gases were considered calorically perfect. This assumption permits closed- 
form solutions while maintaining a physical model reasonably close to that which would be 
observed. Such solutions are invaluable in the assessment of trends and the emphasizing 
of salient points. 

(2) All shock waves with the exception of the adjusting shocks, were assumed to be 
strong. This is a commonly made assumption in shock-tube work resulting in more fac- 
ile solutions with little error. 

(3) Shock wave and interface speeds through nozzles with arbitrary A were taken 
to be invariant. In addition, it was assumed that nonsteady interactions between shock 
waves and trailing entropy discontinuities, due to passage through nozzles, were negli- 
gible. Without these assumptions, analytical solutions of the passage of waves and dis- 
continuities of various types through nozzles would be impossible to obtain. It was felt 
that errors which must result from such assumptions may be held to acceptable levels if 

n = o(lO"l) and jls = o(lQ~l). This results in average wave and inter- 
face velocities relatively unaffected by variations within the nozzle. Trends sought and 
determined under these restrictions are felt to be reasonably accurate. 

(4) Conditions 5' are approximately equal to conditions 5 for configuration m. 
This is an outgrowth of assumption (3) and is necessary for an analytical solution for 

l s 2 /Z n . It was shown that the error which may result from this assumption was not 
enough to invalidate the resulting conclusions. However, as an indication of the extremes 
of this error, figure 23 presents the ratio of the pressure ratios across an unsteady 
expansion followed by a steady expansion to that across the reverse process. In each 
case, the initial and final Mach numbers and A of the nozzle are identical. From this 
figure, it can be concluded that the assumption for preliminary information is justified. 

(5) It has been assumed that the flow in all nozzles is one dimensional. This 
assumption allows analytical solutions throughout the nozzle and is reasonable as long as 
it is not used to describe expansions through large A in extremely short nozzle lengths. 
When such a situation exists, errors in flow properties may become significant. Account 
must also be taken in such cases, of the shock waves necessary to turn the flow through 
the angle that will be formed at the junction of the nozzles and the acceleration chamber. 
These standing waves may interfere with the test flow to the extent that it becomes use- 
less. The configuration producing such a flow would be of no value. 

(6) It has been assumed that the unsteady expansion fan, after passing through the 
nozzle, may be treated as a centered expansion. It is clear from the general unsteady 
characteristic equations (ref. 4) that this is an approximation. It is believed, however, 
that because of the extent of the expansion fan downstream of the nozzle, and with the 
determination of the "virtual center" of expansion, the effects of this approximation can 
be minimized and reasonable results obtained. 
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Configuration I 


From the foregoing analysis, it has been determined that for A = 10 and 
2 ^a^agj , the perfectly expanded configuration suffers severe penalties in 
driver pressure requirements when compared to the basic expansion tube. 


Figure 5 indicates that the driver sound speed requirements for a perfect expansion 
are modest - that is, (a^/agj^- = 11.0 for A = 10, Mg ^ 50. These requirements are 
well within the realm of present-day technology and pose no problems. 


For configurations I and II, the intermediate -chamber parameter, ] for 

p 5\ a l/ 

arbitrary values of A, is represented in figure 13 by the curve for A = 1.0. The ratio 
Pl/Pg, which for this case is identical to that for the expansion tube, increases with Mg 


and decreases with a 


5/ a l- 


The acceleration-chamber parameter 


P 5 f a 10 N 
p 10\5 


for configurations I, n, and III is 


shown in figure 14 and is independent of A. The pressure ratio Pg/P^g increases with 
Mg, decreases with a^g/ag, and, being independent of A, is identical to the same param- 
eter computed for the expansion tube. 


The required driver length for most efficient operation of configuration I is pre- 
sented in figure 15 and is compared at the same values of a^ with its counterpart for 
constant area tube operation, this being denoted by the curve labeled "Expansion tube." 
The comparison reveals that even for relatively long nozzles, that is, l n jl s ^ ~ 0(1), 
advantages may be realized when A > 40. 


The acceleration-chamber length required for maximum test time is shown for con- 
figurations I and II as well as for the expansion tube in figure 16. As previously noted, 

this length is independent of A and is, for reasonable values of , quite large. So 

1 

large in fact is l s , that it must necessarily be made shorter than the length yielding 

Z 

maximum test time even though use is made of only a fraction of the processed gas. Eco- 
nomic as well as viscous interaction problems force this decision. 


Configuration II 

The driver pressure parameter 4 >jj is shown in figure 10. From this figure, it is 
apparent that this configuration requires a driven pressure substantially higher than that 
necessary for the expansion tube. Comparisons are made at equal driver sound speeds. 
There are, however, significant reductions in p^ for configuration II when compared with 
those for configuration I. These gains come at the expense of driver sound speed a^ 
being higher for configuration II. 
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Limitations on the operating range for configuration II are, at the higher end, values 
of Mg at which a perfect expansion occurs and at the lower end, values of M g such 
that M e g M es . These and the possibility of a boundary-layer-separation limitation have 
been discussed. 

The comparison of the driver length of configuration II with those of configuration I 
and the expansion tube is made clear by a consideration of figures 15 and 17. The driver 
for configuration II is considerably longer than that for configuration I and slightly shorter 
than that for the expansion tube. This difference is a result of the higher sound speeds 
required for configuration II and the differing expansion processes in the facilities. 

Intermediate- and acceleration- chamber pressures as well as acceleration- chamber 
length were included in the discussion of configuration I. 


Configuration III 


The driver parameter is shown in figure 12. At the lower values of Mg, 

configuration in permits some reduction of p^ when compared with configurations I 
and H. This, however, is true primarily for higher values of a^. For low driver sound 
speeds, the limiting state for the expanding driver gas is reached at a lower value of Mg 
than that for the expansion tube, resulting in higher driver pressures at moderate values 
of test-section Mach numbers. Configuration in, with regard to p^, is used to best 
advantage with large values of a^/ag as seen from figure 12. 


The intermediate -chamber parameter 
by figure 13. From this figure, the value of 

a 5/ a l* 


1 / 

-2- is portrayed for configuration HI 

^ \3'i 


and M 


5 

would 


P 5 \ 

Pj is seen to increase with A 

and to decrease with the square of ag/aj. One result of such an increase in pj^ 
be the inhibition of boundary-layer effects and hence, the "leaky piston" effect reported in 
reference 5. A consequence of this would be that considerably more flow would be avail- 
able which would be suitable for processing to test conditions. 


The quantity 




is presented in figure 18 and lies within the range of 2.0 


to 0.001, the low end of the range occurring at high Mg and low The param- 

eter is seen to be a strong function of Mg and A. 


Configuration EH now appears to offer the greatest possibilities as far as reducing 
the driver length is concerned. This conclusion, however, will bear close scrutiny since 
for configuration I values of j^s^j ~ o(lO _ 2) can be achieved which are independent 


of Mg, whereas this parameter increases rather rapidly with decreasing Mg for con- 
figuration III. In any event, configuration 1H allows a driver, for all sound speeds, which 
is considerably shorter than that required for the expansion tube; this becomes clear from 
figure 18, the expansion tube curve being denoted by A = 1.0. The reason for this trend 
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is simply that in configuration III, the reflected u + a wave and the interface are allowed 
to intersect earlier (at the end of the driven section) than in the expansion tube. This is 
manifested by a shorter driver. 


The critical component of configuration III then appears to be the acceleration 
chamber. The mass flow analysis described previously is used and the parameter 
(Z S 2 /f n ^ evaluated and shown in figure 20. As pointed out, xjs is the ratio of mass 
unusable for testing purposes to the mass which has flowed past the test section (at the 
end of the expansion or accelerating section) in the interval between the arrival of the 
secondary interface and the last characteristic of the secondary u - a expansion fan. 

K l s ^jl n is chosen such that ip < 1.0, usable flow exists at the test section. With this 
in mind, figure 20 makes it clear that a lengthy acceleration chamber is required if gra- 
dient free test flow is desired. Such excessive lengths may prohibit the use of configura- 
tion III. 

Finally, figure 22 presents the maximum nondimensional nozzle lengths. In config- 
urations I and II the nozzle length parameter l n jl s played a role only to the extent 
that it influenced the driver length parameter For configuration III, (injls^j^ 

is considerably more important in that values of this parameter greater than those shown 
in figure 22 result in zero test time. Consequently, nozzles used must be shorter than 
indicated in figure 22 and this, since (lnjls , ^ j * s ra -ther small in most cases, results 

in a severe violation of the assumption of one-dimensional flow in the nozzles. The most 
severe consequence of this violation, not only for configuration III but also for configura- 
tions I and II, is the formation of recompression shock waves when the flow, which is not 
one dimensional, exiting from the nozzle is collected in a cylindrical section such as the 
intermediate and acceleration chambers. Any such strong shock waves will render config- 
urations in which they occur useless. 


CONCLUDING REMARKS 


The configurations presented would lack feasibility in comparison with the config- 
uration presented in NASA Technical Report R-223. This conclusion is supported pri- 
marily by a consideration of the component lengths, driver pressure, and the deviation of 
the flow in the nozzle from the assumed one-dimensional flow. 

Configuration I is rejected due to the extreme driver pressures required (although 
with relatively low sound speeds). 

Configuration II is rejected as a result of high driver pressure requirements coupled 
with the severely truncated operating regime. 


32 



Configuration in, which was believed to be most feasible, is rejected primarily 
because the short nozzle lengths required in order to insure the availability of testing 
time. Use of such short nozzles would result in strong nozzle recompression shock waves 
disrupting any useful test flow. It should be pointed out that similar standing waves will 
be present to some extent in configurations I and II possibly creating serious flow disturb- 
ances. Shortcomings of this nature are not present in the facility of NASA Technical 
Report R-223 since the test section is located immediately after the expanding nozzle, 
the test flow being utilized prior to its interaction with the recompression waves. Finally, 
it may be argued that configuration in requires higher driver pressures than those in the 
expansion tube (with the same test- section diameter) with little improvement in the viscous 
problems associated with the acceleration chamber. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 4, 1965. 
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APPENDIX A 


ADJUSTING SHOCK ANALYSIS 
The following general flow situation exists (sketch Al): 


<?- 


u l» a l 


u 


es 


u. 


u l >u 2 


Sketch Al 

where u^, Ug, and a^ are known and the shock of Mach number 


M es 


U 1 +u es 


(Al) 


compatible with this situation is desired. From normal shock relations (ref. 6) in shock 
fixed coordinates the expression 


is determined where 


Defining 


gl _ P2 = (y + l)M es 2 

^2 P 1 M es 2 (y - 1) + 2 


u. = u. + u 
1 1 es 


u 2 = u 2 + u es 


? = 


u l-u 2 


(A 2) 


(A3) 


(A4) 
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and substituting this and equations (A3) into equation (A2) gives, after some manipulation, 
the desired shock Mach number as follows: 



APPENDIX B 


NOZZLE LENGTH RATIO FOR CONFIGURATION IH 


The nozzle length which results in zero test time for configuration IH is to be deter- 
mined. Consider first sketch 3 which is a schematic representation of the situation 
described in this appendix. The nozzle length in this case is such that the reflected 
u + a expansion characteristic, the last characteristic of the secondary fan, and the first 
usable test particle coalesce as all reach the end of l S g In order to analyze the 
arrangement of waves shown in sketch 3, a parameter a must be introduced in order to 
permit the investigation of the unsteady expansion exiting from the nozzle, with equations 
derived for a centered u - a expansion in simple flow. This parameter represents the 
distance the first and last characteristics of the expansion must be extended past the 
nozzle exit, in the direction of decreasing x (measured from the nozzle exit), in order 
to intersect or reach a "virtual center of expansion." It is apparent that the center of 
this effective expansion is both an approximation and not at the nozzle entrance since the 
characteristics of the fan have been accelerated through the nozzle. Consequently, a 
is a working parameter, undetermined, and eliminated in the solution of the equations. 

In sketch 3, the ratio of exit to entrance distances for the particle path shown is 


a + l s 2 
a 


M 5 - 1 
M e - 1 



(Bl) 


where equation (68) has been used. For the u + a wave, the ratio of exit to entrance 
distances is given by 


a+ l s 2 

<* + x e 


M 5 - 1 

M e 1 



y - I 

2 

y - 1 

2 


3=Z 



(B2) 


which has been derived by applying equation (36) across the unsteady fan separating 
regions e and 5, multiplying both sides by (ug - a^y^Ug - a e j and substituting equa- 
tion (la) of the present paper and equation (11) of reference 7. 
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Two expressions for r_ 


T e X e\u e -a e u e + a 


(B3) 


T e -^ i + fe (« + -0 


(B4) 


are determined where x e is the distance from the beginning of l (or end of the 
nozzle) to the point at which the u + a characteristic enters the fan and 


T . = JLL r(y - 11 r - 1 

• 1 r» \ ^ 


1 a 


l 2 t 2 y + 1 

The accelerating chamber length l s is given by 


(B5) 


fl B r 


1 s 2 U ^ ln 


(B6) 


the tilde notation indicating a value of xjy equal to unity. 

If the variables in equations (Bl) to (B6) are eliminated and equation (lb) is used 
across regions 2 and e, the following expression for s ^j results: 


-n 


s l> 


(y - l)(y + l) 


2(1 M e 2 


1/2 


y 


I]i(M e + l)(M e - 1) 


m 


s 2 




m 5 - 1 


r l + 


y 


- l 


„ 

2 5 


y + 

y+l 

2(y-l) 


M e “ ^ +L -^ M e/ 


- 1 


M 5 - 1 1 

/ 1+ r^.M 5 \ 

M e - il 

b+Z-LLit] 


/-I 


- 1 


1(«- - «j( M e 2 - >) 


(B7) 
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APPENDIX B 


Equation (B7) is a function of A and M 5 and gives the nozzle length for which zero 

test time exists at the end of • Values of 

b 2 

equation (B7) produce usable test time. 


f l n ll s \ smaller than those given by 
V > VIII 
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Time, t Time 


Shock Wave 



(a) Configuration I. 



(b) Configuration II. 
Figure 1.- Wave schematic. 
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(b) Configuration III. 
Figure 2.- Physical schematic. 





Ratio of driver to test-section pressure, (p*/Pc) 



Figure 4.- Variation of driver pressure ratio with test-section Mach number for basic expansion tube. 
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Asymptote for a./a^ = 5 curve 














o 10 20 30 40 

Test-section Mach number, M 5 


Figure 9.- Adjusting shock Mach number required for configuration 1 1. 


Driver pressure parameter. 



(a) (a 4 /a^_ = 5 and 10. 

Figure 10.- Variation of driver pressure parameter for configuration II. 
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Figure 10.- Concluded. 
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Asymptote 


20 30 40 

Test-section Mach number, 


(a) A = 1 and 10. 

Figure 12.- Variation of driver pressure parameter for configuration 1 1 1. 



Test-section Mach number, M^ 


(b) A = 10 2 . 


Figure 12.- Continued. 



Driver pressure parameter, 



(c) A = 10 3 . 
Figure 12.- Concluded. 
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Pressure parameter, P-j/Pg ( a^/ 



Test-section Mach number, 


Figure 13.- Variation of intermediate-chamber pressure parameter for configuration III. Curve for A - 1.0 is also valid 

for configurations I and 1 1 for all A. 
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Figure 14.- Variation of acceleration-chamber pressure parameter for all configurations. 




Length ratio 



Figure 16.- Ratio of acceleration to intermediate chamber length for configurations l and 1 1 







Test-section Mach number, Mg 

,b) ( a 4/ a 5)* = 10 -°- 


Figure 17.- Continued. 
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Test-section Mach number, M^ 


(b) (a 4 /a 5 j_ = 10.0. 


Figure 18.- Continued. 







Distance, x 


Figure 19.- Wave schematic of nozzle and acceleration chamber for configuration I II. 








Test-section Mach number. Mg 


Figure 22.- Ratio of maximum nozzle length to acceleration-chamber length for configuration 1 1 1 


Ratio of expansion ratios, (p 2 /p 5 ) u _s/(p 2 /p 5 ) 
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